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ABSTRACT 

Scaling relations of metallicity (O/H), star formation rate (SFR), and stellar mass 
(M s tar) give important insight on galaxy evolution. They are obeyed by most galaxies 
in the Local Universe and also at high redshift. In a companion paper, we compiled 
a sample of ~4100 galaxies from redshift to > 3, spanning almost two orders of 
magnitude in metal abundance, a factor of ~ 10 6 in SFR, and of ~ 10 5 in stellar 
mass. We have characterized empirically the star-formation "main sequence" (SFMS) 
and the mass-metallicity relation (MZR) for this sample, and also identified a class 
of low-metallicity starbursts, rare locally but more common in the distant universe. 
These galaxies deviate significantly from the main scaling relations, with high SFR 
and low metal content for a given M sta r- In this paper, we model the scaling relations 
and explain these deviations from them with a set of multi-phase chemical evolution 
models based on the idea that, independently of redshift, initial physical conditions 
in a galaxy's evolutionary history can dictate its location in the scaling relations. Our 
models are able to successfully reproduce the O/H, M sta r, and SFR scaling relations 
up to z £ 3, and also successfully predict the molecular cloud fraction as a function of 
stellar mass. These results suggest that the scaling relations are defined by different 
modes of star formation: an "active" starburst mode, more common at high redshift, 
and a quiescent "passive" mode that is predominant locally and governs the main 
trends. 

Key words: galaxies: abundances - galaxies: dwarf - galaxies: evolution - galaxies: 
high-rcdshift - galaxies: starburst - galaxies: star formation 



1 INTRODUCTION 

Scaling relations of metallicity O/H, stellar mass M s tar, 
and star-formation rate (SFR) are powerful tools for 
probing galaxy evolution across cosmic time. The star- 
formation "main sequence" (SFMS) co nnects current star 
formation with a g alaxy's stellar mass (iNoeske et al.ll2007l: 
Salim et all 120071 ; ISchiminovich et all 120071 ; iBauer et all 
20111 ). and the mass-metallicit y relation (MZR) link s stel- 
lar ma ss with metal c o ntent |Tremonti et al.l |2004 ). Re- 
cently, IMannucci et al.l (|2010h found that SFR influences 
a galaxy's position in the MZR, and introduced the "Fun- 
damental Metallicity Relation" (FMR), which reduces the 
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scat ter in their local samp le of galaxies to ~0.06 dex (see 
also lLara-Lopez et all |2010T ) . 

There is, however, evidence that galaxies outside the Lo- 
cal Universe (z > 1) do not follow the same SFMS and MZR. 
The FMR successfully fits Lyman Break Galaxies (LBGs) 
up to z < 2, but fails for higher redshifts (|Mannucci et al.l 
l201Cf ). The SFMS has similar slopes w ith increasing red- 
shift but larger offsets (|Bauer et al.ll201lh ; for a given stellar 
mass, a more distant galaxy has a higher SFR. Starbursts at 
all redshifts tend to lie above the SFMS, with high specific 
SFRs (sSFRs), and to o much luminosity (or stellar mass) for 



their metallicity (e.g., Hovos et al.ll2005l : iRosario et al ] |2008l : 
ISalzer et aI.ll2009l ; |Peeples et alJl200^ TThese high-SFR out- 
liers have been attributed to merger-induced episodes of 
star formation that tend to be short-lived and intense (e.g., 
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iRodighiero et al.ll20lil ). while the SFMS is followed by more 
passively evolving galaxies that form stars over longer times. 

Here we propose that mergers may not be directly re- 
sponsible for producing the starbursts that depart from 
the general SFR-mass-metallicity relations. Instead, inde- 
pendently of their cause, different modes of star forma- 
tion arising from different initial sizes and densities, may 
be driving the deviations. Active star formation that takes 
place in compact regions with dense gas has distinct prop- 
erties from a more passive mode in larger, more tenuous 
regions; in the former, dynamical times are shorter, infrared 
luminosities and S FRs are higher, and mol ecular fractions 
are more elevated (|Hirashita fc Hunt] [200jj) . In the Local 
Universe, compact sizes and high-density gas characterize 
metal-poor starbursts which occur in dwarf galaxies with 
high SFRs (|Hunt fc Hirashitall2009T ). At high redshift, star- 
bursts also tend to be compact in size (< 4 kpc in diameter, 
iTacconi et ail 120061 . l20Qg| ; iToft et all |2009| 1H Compactness 
also shapes spectral energy distributions causing them to 
peak at shorter wavelengths, because t he dust heated by a 
compact starburst tends to be warm er jChanial et al.ll2007l ; 
iGroves et alj|2008l ; [Melbourne et aT1l2009l ). 

Active star formation is expected to also be charac- 
terized by high ionization parameters because of the com- 
pact and dense nature of the active SF complexes. Fig- 
ure [T] shows the ionization parameters U of a sample of 
low- metallicity Blue Compact Dwarf galaxies (BCDs) taken 
from [Hunt fc Hirashital (|2009l '). The ionization parameter of 
the Hn regions increases with ionized gas density n e ; for this 
sample the significance of the correlation is >99%. Because 
the s ize of the Hn region is in versely correlated with density 
(e.g-. lHunt fc Hirash ita 2009), U depends almost linearly on 
n e . Anomalous Hn-region excitation is freq uently found in 



some high-redshift gal a xy populations (e.g ., Hammer et al 
1997I ; iLiu et all |200St iBrinchmann et alj |200g| ; lErb et a" 



20101 '). Excluding a contribution from a weak active galactic 



nucleus, it is likely that this excitation results from the ex- 
treme physical conditions of the active SF mode, with dens e 
ionized gas and high SFR surface densities |Liu et alj|200ct ). 
Such properties naturally arise in dense and compact Hn 
regions with particularly hard radiation fields. Thus, the 
anomalous excitation and high ionization parameters ob- 
served in high-redshift samples may also be signatures of 
active starbursts which occur in a chemically unenriched in- 
terstellar medium (ISM). 

Besides compact, dense star-forming regions, and possi- 
bly high-ionization parameters for the nebular gas, there are 
further, probably related, properties shared by some galaxy 
populations in the n earby universe and those at high red- 
shift. As discussed in lHunt et al.l |2012l . hereafter Paper I), 
some metal-poor BCDs deviate significantly from the SFMS 
and the MZR; they are characterized by high SFRs, high 
sSFRs, and are situated in the SFMS and the MZR in the 
same region as the z — 3 LBGs. For a given (sub-solar) 
metallicity, these galaxies can have an excess of stellar mass 
of two orders of magnitude, and a factor of ~ 10 excess 



1 However, quiescent galaxies at z ~ 2 tend to be more compact 
than starbursts as the same redshift, implying that those galaxies 
were formed by compact nuclear starbursts at even higher red- 
shifts, 2 ~ 3 — 4. 
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Figure 1. Ionization parameter vs. root-me an-square electron 
densit y of Hn regions in a sample of BCDs from lHunt fc Hirashital 
teOOgh . The typical error bars are shown in the lower right corner. 
The dotted line shows the best-fit correlation. 



SFR for a given mass. Some of the "Green Peas" and Lumi- 
nous Compact Ga l axies (LC Gs) at z ~ 0.1 — 0.4 selected by 
ICardamone et all (|2009l ) and llzotov et alj (|201lh also lie in 
the same regions of the SFMS and the MZR and show simi- 
lar excesses. Similar physical conditions seem to occur at all 
redshifts, but different galaxy populations are observed at 
different redshifts because of how they are selected. Since ex- 
treme starbursts are rare locally, but more common at high 
redshift, it could be that different modes of star formation 
drive the scaling relations of stellar mass, metallicity, and 
SFR we observe. 

In this paper, we present semi-analytical models which 
are aimed at understanding how different modes of star 
formation affect the SFMS and MZR. In particular, we 
focus on metal-poor galaxies with high SFRs, which are 
generally more massive than would be expected given 
their metal abundance. We develop three families of ini- 
tial conditions for star formation and compare model pre- 
dictions with several samples comprising < 1100 nearby 
metal-poor dwarf galaxies, intermediate redshift and high- 
redshift galaxy populations (see Paper I). The initial 
conditions are designed to embrace the active/passive 
(compact+dense/extended+diffuse) modes of star forma- 
tion and thus distinguish the active starburst modes from 
more passive quiescent ones. The questions we aim to ad- 
dress are: 

(a) Why can both high and low SFRs be observed at the 
same metallicity? 

(b) What are the specific conditions which distinguish star- 
bursts (active SF mode) from quiescent star formation 
(passive mode)? 

(c) Is there a connection between active dwarfs and high- 
z starbursts in terms of their scaling relations (SFMS, 
MZR, FMR)? 

The paper is organized as follows: in Sect. [2] we briefly de- 
scribe the samples we have compiled and how we calculate 
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the stellar mass. The sample is more completely described 
in Paper I. In Sect. [3] we discuss our semi-analytical mod- 
els which consider different sets of initial conditions in order 
to mimic the starburst and quiescent modes of star forma- 
tion. Comparisons of model predictions with observations 
are given in Sect. 2] and Sect. [S] discusses the relation of the 
gas-to-stellar mass ratio and metal abundance. In Sect. [6] 
we describe possible similarities between nearby and distant 
starbursts, and how the scaling relations we observe and 
successfully model can shed light on how galaxies evolve. 



2 THE SAMPLE AND THE OBSERVABLES 



Our sample was compiled on the basis of three "pseudo- 
observational" parameter^]: SFR, metal abundance [as de- 
fined by the nebular oxygen abundance, 12+ log(0/H)], and 
stellar mass, M s tar- We included only galaxies with these 
data either already published, or for which we could derive 
the quantities ourselves (in particular the stellar masses). 
These criteria were met by galaxies from five samples at 
z ~ 0: 21 ne a rby d warf irregular galaxies (dlrr) published 
by iLee et all (|200rl hereafter "dlrr"), 129 of the HMpc 
distance-limited samp le of nearby galaxies (11HUGS, LVL: 
iKennicutt et al.| 20081). and 89 B C Ds in the samp l es pre - 
sented bvlEngelbracht et all (|2008l ), lFumagalli et all (|2010t ). 
and lHunt et al.l (|2010l ). Seven galaxy samples at higher red- 
shifts also have the necessary data: the "Green Pea" compact 
alaxies identified by th e Galaxy Zoo team at z ~ 0.1 — 0.3 
Cardamone et ai1l2009l). the LCG sample at z ~ 0.1—0.6 se- 



lected by llzotov et al.l d201lf). LBGs a t z ~ 1 JShaplev et al.1 



)ley e 

2 (IShapIev et alj|2004 lErb et all 12009 ). and 



12005 j ). z 

z ~ 3 (|Maiolino et"~afl 120081 ; rMannucci et~afl 120091 ). When 
stellar masses were not available in the literature, we deter- 
mined stellar masses from IRAC photometry as discussed 
below, and in detail in Paper I. 

Altogether, we consider in the analysis 1070 galaxies 
from z ~ to z ~ 3 with the requisite three quantities 
of SFR, 12+ l og(Q/H), and M Btar . Excep t for the galaxies 
with z > 1 jShaplev et al.1 12004|. l2005al; lErb et al.l 120061 ; 



iMaiolino et al l |200S|; iMannucci et al. 1 12009| ~ and a few of 
the LVL galaxies (jMarble et al.l 120101 ) . metallicities for all 
galaxies are deter mined by the "dire ct" (electron tempera- 
ture) method (e.g.. Ilzotov et al]l2007l ) since the samples are 
defined by requiring detections of OinA4363. SFRs are de- 
rived from Ha luminosities corrected for extinction; in a few 
BCDs, the IR luminosity gave higher SFRs, so we used the 
higher value. SFRs using the iKennicuttl (1 19981 ) conversion 
are reported to the IChabrierl ( 20031 ) Initial Mass Function 
(IMF). Paper I gives a detailed description of the sample. 

We calculated the stellar masses for the local samples 
(dlrr, BCDs, 11HUGS/LVL) from Spitzer/IRAC observa- 
tions at 4.5 um (and when a vailable al s o 3.6 um). We based 
our method on that used bv lLee et al.l J2006I) who analy zed 
sub-solar metallicity models from Bell fc de Jond (|200ll ) to 
derive stellar mass-to- light (M/L) as a linear function of 
B — K. The advantage of using 4.5 urn to measure stellar 

2 We say "pseudo-observational" because these variables are de- 
rived from observations rather than being directly observed. Nev- 
ertheless, for simplicity, we will call O/H, SFR, and M a t a r "ob- 
scrvables" . 



masses is that it minimizes variations in the M/L ratio (e.g., 
lJun fc Imlliiosl ). because the ratio depends only weakly on 
age and metallicity, and dust extinction is negligible. 

However, before calculating the masses, we subtracted 
the nebular continuum and recombination line emission 
from the broadband 4.5 vim photometry. Nebular emission 
in metal-poor starbursts, both in the co ntinuum and lines , 
can significantly affect photometry (e.g.. iReines et alj|20ld : 
lAtek et alj|201ll) . and thus have a potentially strong impact 
on the inferred stellar mas ses. This can be especial ly impor- 
tant in the 4.5 um band dSmith fc Hancockfl2009l ). Hence, 
before applying the formalism of lLee et al. [ ~(|200aV from the 
SFR we have inferred the strength of Bra emission and 
the nebular continuum in the 4.5 urn IRAC band , and sub- 
tracte d it from the observed flux. Then, following ILee et all 
|2006l ). 4.5 [im luminosity was converted to stellar mass with 
a M/L corrected for different ages and metallicities with a 
color correction B— [4.5]. All stellar masses are scaled to the 
IChabrierl |2003l ) IMF. 

The final combined sample covers a ~ 10 J range in 
stellar mass, a factor of ~ 10 6 in SFR, and 2 orders of 
magnitude in oxygen abundance. The metallicities are gen- 
erally (formally) quite accurate, with uncertainties on the 
direct method of 0.05 dex, but the strong-line method 
and its comparison with the direct method of determining 
abundances from electron temperatures are compromised 
by differences in calibr ation that can be quite large (e.g., 
iKewlev fc Ellisonll2008l . see also Paper I). The SFRs suffer 
from the typical uncertainties of Ha measurements and the 
uncertain extinction corrections, and the uncertainty is ex- 
pected to be £ 30%. The uncertainties on the stellar masses 
are typically 0.1 dex (30%), but can be as high as 50-70% 
when there are unusual colors. In any case, the large dy- 
namic range of our combined sample should overwhelm any 
uncertainties on individual galaxy parameters. More details 
are given in Paper I. 



2.1 Metal-poor starbursts 

Since we are interested in deviations from the scaling re- 
lations, particularly those of galaxies with high SFR, rela- 
tively high mass, and low metallicity, in Paper I we quanti- 
fied a "low- metallicity starburst". Because deviations from 
the MZR become significant at abundances at or below 
12+ log(O/H)~8.0, we consider this as a threshold below 
which a galaxy can be considered "low- metallicity" . De- 
viations from the SFMS and the MZR become notice- 
able at SFRs of ~ 0.6MQyr~ 1 . Hence we loosely define 
low- metallicity starburst (LMS) as a galaxy with a SFR 
>0.6M o yr _1 aixd with 12+log(O/H)<8.0. This is a rather 
arbitrary label, but one that enables us to identify the galax- 
ies most likely to deviate from the scaling relations (see Pa- 
per I). Among the various samples, the LMSs are relatively 
rare: 2% of the LVL fall into this category; ~11% of the 
BCDs; 22% of the LCGs; and ~27% of the high-z LBGs. 
Clearly selection effects (e.g., emission line surveys for the 
BCDs, emission-line flux for the LCGs) make it more likely 
for a galaxy to be a LMS, but there is some indication that 
LMSs are more common at high redshift. 
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3 THEORETICAL MODELS 

Here we develop the models with which we explore the idea 
that particular initial physical conditions in a galaxy's evo- 
lutionary history can change its location in the scaling re- 
lations independently of redshift. In particular, we have de- 
vised three sets of starting parameters which dictate the sub- 
sequent evolution of the galaxy, and which encompass the 
observed distinctions between "active" and "passive" star 
formation (see Sect.QJ. The models and the different initial 
conditions are described below. 

Although our models have been shown to realistically 
approximate the chemical evolution in nearby galaxies (e.g., 
Molla et all 1 19961 . Il997l : IMolla fc Di'azl 120051 : iMagrini et ail 



20071 ). in those models gas accretion has been considered. 
Here, as a first approximation, we use a closed box in which 
neither infall of pristine gas nor outflow of metal-enriched 
material are contemplated. Despite these significant limita- 
tions, which will be addressed in future work, we will show 
below that our models are remarkably successful in predict- 
ing the details of the SFMS and MZR scaling relations. They 
also are consistent with the observed gas fractions both lo- 
cally and at high redshift. 



3.1 The multiphase chemical evolution model 

The models developed for the present work have been re- 
named as MICE (Multiphase Interstellar medium Chemical 
Evolution) models. MICE models ado pted here are a gener - 
alization of the multi- phase model bv lFerrini et al.l (|l992h . 
originally built for the solar neighb orhood, and subse - 
quently extended to the entire Ga laxy jFerrini et al.lll994[) . 
and to other disk galaxies (e.g.. IMolla. et all 1 1999. Il997l: 



IMolla fc Di'azll2005l : lMagrini et al.l l2007. 200<1). We refer to 



those papers for a detailed description of the general for- 
malism. The main differences of the models presented here 
relative to those used to describe the chemical evolution of 
disk galaxies are: i) MICE models are one-zone models, with 
a single galactic component and without infall or outflow 
processes; ii) no radial variations are considered, so that a 
uniform evolution is assumed for the whole galaxy. 

Each galaxy is composed of four components: diffuse 
gas (g), clouds (c), stars (s) and stellar remnants (r). At 
t = all the mass of the galaxy is in the form of diffuse 
gas; subsequently, the mass fractions are modified by sev- 
eral conversion processes. The main processes considered in 
the model are: (i) conversion of diffuse gas into clouds; (ii) 
formation of stars from clouds and (Hi) disruption of clouds 
by previous generations of massive stars; (iv) evolution of 
stars into remnants and return of a fraction of their mass to 
the diffuse gas. 

We describe in the following the main equations of the 
model. Stars form by cloud-cloud collisions with a rate P s f 
and by the interactions of massive stars with clouds with 
rate Pi (induced star formation), 



±l± = P sf f 2 c+ P t f s f c ~Df a 



(1) 



where f s (t) and f c (t) are the mass fraction in stars and 
clouds, respectively, and D is the stellar death rate. The 
quadratic dependence of the SFR on the mass of clouds re- 
sults from the coUisional nature of the process. Clouds con- 



dense out of diffuse gas with a rate P c and are disrupted by 
cloud-cloud collisions and by winds from massive stars, with 
rates Pd and P w , respectively, 



p f3/2 



{p S f + Pd)f;-(Pw + Pi)fsfc 



(2) 



die 
dt 

where f a (t) is the mass fractions of diffuse gas. The 3/2 
exponent results from the assumption that the process of 
condensation of clouds occurs on the time scale of the grav- 
itational instability in the diffuse gas, proportional to the 
inverse square root of the gas density. This introduces an 
extra power of 1/2 in the gas mass, because mass fractions 
are equivalent to mass densities given our assumption of 
constant galactic volume (see below). The mass fraction of 
diffuse gas evolves as 



d/ 9 _ D r 3/2 , 



dt 



= -Pcft* +Pdfc+P W fsf. 



(■■>>) 



where we have assumed that clouds disrupted by cloud-cloud 
collisions and stellar winds produce diffuse gas. Finally, stel- 
lar remnants are produced according to 

dfr 



dt 



= Df s 



(4) 



Thus, the total mass of the galaxy does not change with 
time (f g + / c + f s + f r = 1, closed-box model). 



3.2 Model parameters 

At time t = galaxies are represented as uniform-density 
spheres of neutral cold gas with mass density po from which 
we obtain the initial radius of the galaxy 



i? gal = [3Af gal /(47r /5o )] 1/3 . 



(5) 



In the MICE models, galaxies evolve at constant volume, i.e., 
they maintain the initial radius given by Eq. ([5]). In particu- 
lar, we consider two values of the initial densities of atomic 



gas, po = 2 x 10 24 g cm 3 and po 



1 x 10 22 g cm 3 (cor- 



responding to numerical densities of 1 cm -3 and 400 cm -3 ) 
representing slowly-evolving (SE) and rapidly-evolving (RE) 
galaxies, respectively. 

We produced a grid of 30 galactic models by varying 
the following two parameters: 

a) the initial mass: M gaJ = 10 7 , 10 8 , 10 9 , 10 10 , 10 11 M Q ; 

b) the gas density in molecular clouds: umc = 10 3 cm -3 
("diffuse" clouds), 7imc = 5 x 10 4 cm -3 ("compact" 
clouds), and nivic = 10 6 cm -3 ("hyperdense" clouds). 

The molecular cloud (MC) densities are chosen to 
encompass the range of volume de n sities obser v ed in 
MC complexes (e.g., ILarsonl Il98ll : Ijiiina et al. I Il999l : 
IKrumholz fc Tar] 120071 ). The relatively high densities cor- 
respond to the dense cloud cores revealed by high-density 
gas t racers, since these are the true sites of star formation 
(e.g.. IL"ada et al.lll99ll ). Obviously, we cannot capture the 
complexity of star formation with our models, but these 
densities should be representative of the r ange of conditions 
in quiescent galaxies and star bursts (e.g.. lAalto et al1ll995l : 
IKrumholz fc ThompsonhoOTl '). Galaxy mass is varied to as- 
sess whether dwarf galaxies, with similar atomic and molec- 
ular cloud properties as their more massive counterparts, 
can be considered as scaled-down versions of giant galaxies, 



or -equivalently- whether massive galaxies can be viewed as 
scaled-up versions of dwarfs. 

We associate the diffuse clouds with the passive mode of 
star formation, and the hyperdense clouds with the active 
one. The distinction between active and passive modes is 
relatively arbitrary, so the compact clouds defined here oc- 
cupy a "grey area", forming stars in an intermediate mode 
between "active" and "passive". Obviously, star formation 
in real galaxies encompasses a wide range of densities; these 
three classes are a simplification intended to illustrate the 
effects of different physical conditions which could prevail in 
different galaxy populations. 

The RE galaxies in our models correspond to the star- 
burst mode of star formation thought to be responsible for 
elliptical-gala xy assembly at z ^ 1 — 2 (|De Lucia et al.ll2006l ; 
lRenzinill2006l ). On the other hand, the SE galaxies in our 
models are generally imilar to Local Universe spiral and ir- 
regular galaxies. In what follows, we will designate the two 
classes of models by RE galaxies and SE galaxies; these are 
very loosely associated with proto-ellipticals at high redshift 
and spiral disks in the Local Universe. Our intention is to 
distinguish between the two ways of scaling the models, in 
the sense of spheroid proto-galaxies which evolve rapidly and 
disk-like proto-galaxies which evolve more slowly. 

3.3 Scaling the models 

The coefficients P c , P s f, etc. in the evolution equations 
Eq. {TJ)— (J4j) represent "rates" for the outcomes of the various 
processes. Their dependence on the fundamental parameters 
of the models can be estimated from educated guesses on the 
dominant physical processes responsible for the conversion 
of gas in clouds, stars and vice versa, as outlined below. 
However, their specific values cannot be obtained from first 
principles only, but are determined in general from the ap- 
plication of the model to a large sample of data. 

The models adopted here have been originally devel- 
oped to study the chemical evolution of the Galaxy and 
nearby galaxies (e.g iFerrini et al]|l992l . Il994l ; iMolla fc Dfazl 
120051 ; IMagrini et al.ll2007l . 120091 ). Thus, the coefficients reg- 
ulating the various processes have been calibrated to repro- 
duce the wealth of observational data available for the Milky 
Way and Local Group galaxies. In order to apply our models 
to the galaxies of our sample, it is necessary to adopt sim- 
ple yet robust scaling relations for the process probabilities 
in terms of the fundamental parameters defined in Sect 13. 21 
namely the initial mass M ga i and the density of molecular 
clouds riMC- 

We start from the s et of P c and P a f valu es adopted to 
reproduce our Galaxy in IMagrini et ail (|2009l ). For the most 
massive galaxy of our set (10 11 Mq) of SE galaxies, we adopt 
the P c and P s ; values of the MW, while for the RE galaxies 
of the same mass we increase the P c value by a factor 20. The 
factor 20 arises from the higher density of the gas from which 
the proto-galaxy is formed, resulting in a faster conversion of 
atomic gas into clouds and stars, proportional to the inverse 

1/2 

of the free-fall time, i.e., tx n HI . For both classes of galaxies, 
in case of diffuse molecular clouds we normalize the value of 
P s f to the efficiency value adopted for the MW. We thus 
obtain the values of P c and P a f values for the lower mass 
galaxies and for galaxies hosting denser molecular clouds 
with the relations described below. 
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In add ition, we must take into account the well known 
resul t (e.g.. lMatteuccll994l ; rCalura et aHl2009l ; |pipino et all 
I2009T ) that the star formation efficiency increases with galaxy 
galactic mass M ga i. In our models, the SFR efficiency is the 
inverse of th e depletion timesca le of the cloud mass frac- 
tion (see also lCalura et al"1l2009l ). This mass dependence is 
implied, for example, by t he increasing tren d of [Mg/Fe] 
with Mgai in ellipticals fe.g.. lNelan et aLll2005l ) and from the 
faster evolution of la rge disks with respect to smaller ones 
jBoissier et al.ll200ll ). This relation has taken the name of 
chemical downsizing |Matteuccilll994l ). It is also one of the 
conc urring arguments to explain the mass-metallicity rela- 
tion |Calura et all 120091 ). ICalura fc Matteuccil (|2006T ) sug- 
gested that the main features of local galaxies of different 
morphological types may be reproduced by interpreting the 
Hubble sequence as a sequence of decreasing star formation 
efficiency from early typ es to late types, i.e. from ellipti- 
cals to irregulars (see also lGalli fc Ferrinill 19891 . for an early 
interpretation of the Hubble sequence in terms of SFRs). 

In our models we have scaled the rate of cloud con- 

1/3 

densation as P c tx M 'j and the rate of star formation as 

P s /oc M^'n^Q. The scaling with mass corresponds basi- 
cally to a scaling with size, and, at least for the parameter 
Psf, follows naturally from the interpretation of this quantity 
as a collision frequency between clouds. The latter is given 
by the product of the number of clouds per unit volume 
(roughly proportional to M ga i iZZj) times the cloud's cross 
section (assumed the same in all models) and the cloud- 
to-cloud velocity dispersion (proportional to R~^/ 2 ). 
However, the SFR efficiency also depends on the dynamical 
(free-fall) time of the molecular clouds (tx n^ c ), which re- 
sults in values of P s f larger by a factor of ~ 7 and ~ 32 in the 
case of compact and hyperdense clouds, respectively, with 
respect to diffuse clouds. Observing that what we call "dif- 
fuse" clo uds are typical of the MW and nearby quiescent dis k 
galaxies jSimon et al.ll200ll ; [Krumholz fc Thompson! 120071 ') . 
we can just scale the value of P s f adopted for the MW (taken 
as a reference) to obtain the P s f parameter in the other two 
cases. This scaling approximately reproduces the observed 
fraction of gas (~1%) th at forms stars during a free-fall time 
|Krumholz fc Tanl [20071 ). 

In summary, our models consist of two classes of galax- 
ies: slowly- and rapidly-evolving galaxies. For SE galaxies, 
we normalized the P c coefficient of the most massive ones 
(M gaJ = 10 11 M fl ) t o the coefficient of the MW model of 
IMagrini et al.l (|2009l ); for RE galaxies, we adopted the same 
coefficient multiplied by a factor of 20. For both classes of 
galaxies, in the case of diffuse molecular clouds we normal- 
ized the P s f to the value adopted for the MW. We then de- 
rived the P c and P s f coefficients for the lower mass galaxies 
and for galaxies hosting denser molecular clouds by adopting 
the scaling P s /oc M^n^, and P c cx . We emphasize, 
however, that our scalings are largely empirical, and must be 
considered as a rough parametric representation of physical 
processes not included in detail in our model (internal dy- 
namics, mass loss and accretion, etc.). In Table[T]we present 
a summary of our model scaling factors with respect to the 
MW values (P c ,mw and P s /,mw) for each type of galaxy. 
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Table 1. Model scaling factors with respect to the MW values. 



Mgal Pc/Pc,MW Psf/Psf,MW P s f / P s f ,MW P s f/Psf,MW 



(Mq) Diffuse Compact Hyper-dense 



Rapidly-evolving galaxies (njji = 400 cm 3 ) 

10 11 20. 1 7 32 

10 10 9.6 0.46 3.2 15 

10 9 4.4 0.22 1.5 7 

10 s 2.0 0.10 0.7 3.2 

10 7 0.8 0.02 0.1 0.6 



Slowly-evolving galaxies (nui = 1 cm 3 ) 

10 11 1.1 7 32 

10 10 0.46 0.46 3.2 15 

10 9 0.22 0.22 1.5 7 

10 8 0.10 0.10 0.7 3.2 

10 7 0.02 0.02 0.1 0.6 



range 8 iviq <. m <. u m® 60 ivip, <. m <. iuu ivip, , 
which are not included in tables of lChiefn fc Limongil (|2004r ) 
by linear extrapolation of the yields in the mass range 
13 M Q < M < 35 Mfl. The SNIa yield s are taken from 
the model CDD1 bv llwamoto et all |l999f ). 

The choice of different se ts of stellar yie l ds aff ects the 
final results as highlighted bv lRomano et all {{201(f) . In the 
present work, we are mainly interested on the evolution of 
the oxygen abundance, whose nucleosynthesis in stars is rel- 
atively well understood. Thus, with different prescriptions 
for the O production we do not obtain dramatically differ- 
ent results in terms of O/H, although, as demonstrated by 
I Romano et al.l (|20TOt ). the [O/Fe] ratio may vary 0.2 dex at 
solar metallicity, and even more at lower metallicities. 

Another fundamental ingredient is the IMF. Several 
parametrizations of the IMF have been widely used in 
the literature . In t he present version of the code we use 
the IChabriel (|2003h mass function in the mass range 0.1- 
100 M . 



4 COMPARISON OF MODEL PREDICTIONS 
WITH OBSERVATIONS 

We can now address the active and passive SF modes of 
galaxies with similar masses and metallicities. We first illus- 
trate the time evolution of the stellar mass Af star , the SFR, 
and the metallicity 12+ log(0/H). We then discuss the mod- 
els for the observed scaling relations described in Paper I, 
namely the SFMS, or correlation between M sta r and sSFR, 
and the MZR, or mass-metallicity relation. 

4.1 The time evolution of O/H, SFR, and AW 

The time evolution of metallicity, SFR and M atar are shown 
in Fig. [2] for two sets of galactic models with different ini- 
tial masses M ga i = 10 8 (shown as dark violet/green, for 
SE/RE, respectively) and 10 11 M (blue/red, SE/RE). For 
each mass, we show MCs with three different densities: dif- 
fuse (shown as dotted lines), compact (dashed lines), and 
hyperdense (solid lines). 

For a given galactic mass, M ga i, RE galaxies evolve 
much more rapidly than SE galaxies. They accumulate met- 
als and build a considerable fraction of their stellar mas^f] in 
a shorter time (~ 10 8 yr for Af ga i = 10 11 M Q , and ~ 10 9 yr 
for M ga i = 10 s M©). Slowly-evolving disk-like galaxies of the 
same mass need much more time for metal production and 
to assemble stellar mass (~ 10 9 yr and more than a Hub- 
ble time to produce the same stellar fraction, respectively). 
This signature of chemical downsizing, here closely related 
to mass downsizing, stems from the considerable fraction 
of the stellar mass in massive galaxies formed in less than 
1 Gyr, while 50% or more of the stellar mass of less massive 
galaxies is still being created at the present time. This is be- 
cause of the high efficiency of star formation of RE galaxies 
in their formation phase. 

For a given galactic mass, SFR also peaks at earlier 
times for RE galaxies with respect to SE galaxies, as already 
found by ISandagel (|l986h . Galaxies with denser molecular 
clouds (i.e., a more active SF mode) have a higher maxi- 
mum value of the SFR with an earlier onset than galaxies 
with less dense clouds. Within either SE or RE galaxies, 
for a specific cloud density, more massive galaxies reach a 
higher SFR and at earlier times than less massive galaxies 
(because of their i ncreased star-formation efficiency, see also 
iTosi fc DiaJll985h . 

At a given molecular cloud gas density, jimc, more mas- 
sive galaxies evolve more rapidly than less massive ones 
(e.g., downsizing). For the most massive galaxies, M ga i = 
10 11 Mq, the final oxygen abundances are similar for the 
rapidly- and SE families, with slightly lower values for SE 
galaxies, even though these values are reached over differ- 
ent timescales. For the least massive galaxies with M ga i = 
10 8 Mq, SE and RE galaxies achieve different metallicities at 
the present epoch. Again, this is because of chemical down- 
sizing, since less massive galaxies form stars less efficiently 
Thus, although SE disk-like galaxies with Af ga i = 1O 8 M 
have the potential to reach the same metallicity of RE galax- 
ies of the same mass, they will do so only on timescales much 

3 The stellar mass plotted here is the "luminous mass" , and does 
not take into account the mass of the stellar remnants, white 
dwarfs and neutron stars. 



3.4 Stellar yields and IMF 

The chemical enrichment of the gas is modeled u si ng th e 
matrix formalism developed by iTalbot fc Arnettl (|l973f ). 
The elements of the restitution matrices Qi,j(M, Z) are 
defined as the fraction of the mass of an element j ini- 
tially present in a star of mass M and metallicity Z that 
it is converted into an element i and ejected. This ver- 
sion of the MICE model takes into account two different 
metallicities, Z = 0.02 and Z = 0.006, and 22 stellar 
mass bins (21 for Z = 0.006) from M min = 0.8 M Q to 
M max = 100 Mq, for a total of 43 restitution matrices. For 
low- and intermediate-mass s tars (M < 8Mq) we use the 
yields bv lGavilan et al.l (120051 ) for both values of the metal- 
licity. For stars in t he mass range 13 Mq < M < 35 Mq we 
adopt the yields bv lChieffi fc Limongil (|2004T ) for Z = 0.006 
and Z — 0.02. We estimate the yields of stars in the mass 
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Time [yr] Time [yr] Time [yr] 

Figure 2. Predictions of O/H (left panel), SFR (middle panel), and M s t ar (right panel) vs. time in the MICE models. Models of two 
representative galaxies, 10 s Mq and 10 11 Mq, are coded by line type and color. The diffuse, compact, and hyper-dense cases arc shown 
as dotted, dashed, and solid lines, respectively. Different galaxy masses are distinguished by color: green/dark violet for 10 s Mq, and 
blue/red for 10 11 Mq (SE/RE, respectively). In the middle panel, dotted horizontal lines show the four SFR regimes considered here: 
SFR ^ 0.06 M yr-\ 0.06 < SFR < O.6M yr-\ 0.6 < SFR < 6Mq yr -1 , and SFR ^ eMoyr" 1 . 



longer than a Hubble time. They also require more than a 
Hubble time to fully convert all their gas into stars. 

Within each class (SE, RE), galaxies with the densest 
clouds evolve faster. These differences again result from the 
star formation efficiency which, in our models, depends also 
on the free-fall time of the molecular clouds. 

The left panel of Fig. [2] shows the degeneracy between 
the track of the SE galaxy with M ga i = 10 11 Mq and that 
of the RE galaxy with M ga i = 10 s Mq. The scaling factors 
described in Sect. 13.31 conspire to make the behaviour of 
these two galaxies virtually identical in terms of metallicity 
evolution. However, the evolutionary timescales of the SFR 
and the stellar mass differ, suggesting that metallicity alone 
is not a good indicator of the evolutionary status of a galaxy. 
Specifically, a slowly evolving massive galaxy and a low-mass 
rapidly evolving galaxy will have similar metallicities. 

Indeed, Fig.[2]also shows that for a given galactic mass, 
it is possible to achieve the same metal abundance (left 
panel) with different SFRs (middle) and different stellar 
masses (right). The galaxy can be less massive but, because 
of its active mode, more rapidly form stars, and thus accu- 
mulate more metals and stellar mass on shorter timescales. 
Alternatively it can be more massive and thus produce more 
of everything, but can also form stars and accumulate metals 
more slowly. The key aspect of our modeling is the evolu- 
tionary timescale which is parametrized either through the 
global properties of the galaxies, governed by how we have 
scaled our models (i.e., slow vs. rapid evolution), or through 
the density of the gas in the molecular clouds within the 
galaxies. Either way, shorter dynamical times drive more 
rapid evolution. We will discuss below how these effects can 
help explain the active vs. passive dichotomy, and conse- 
quently how we can better understand the behavior of low- 
metallicity starbursts, and their departure from the scaling 
relations. 

In any case, the present-day abundances of our model 
galaxies are somewhat high. This is almost certainly due 
to the closed-box nature of our models and the lack of 
gas exchange with the environment. In fact, this is what 
led chemical evolution modelers to introduce both infall of 
metal poor gas and later also galactic winds powered by su- 



pernova explosions in order to reduce the predicted metal 
abundances an d bring them to a better agreement with the 
observed ones (|Matteucci fc Tosilll985l ; lTosilll988l ; IPilvuginl 
Il993h . Without galactic winds the products of stellar nucle- 
osynthesis are always retained, and without infall of pristine 
gas they are never diluted. It is well known that gas accre- 
tion and expulsion can sig nificantly alter the chem i cal abun- 
dances in galaxies (e.g., iMatteucci fc Tosil 1 19851 ; IPilvuginl 
1 19931 : iMarconi et"aT]|l994r . iRomano et al.ll200rj ). 

However, infall and/or outflow are not the only mech- 
anisms that could vary the abundances. The absolute level 
of the abundance depends also on the combination of the 
adopted stellar yields and IMF. iMagrini et al.l (|201fj| ) have 
shown that the maximum variation of the abundance level in 
the disk of the spiral galaxy M33 obtained by varying the 
parametrization of the IMF (and keeping fixed the stellar 
yields) is ~0.6dex, with the Chabrier IMF having the higher 
abundances. The impact of adopting different stellar yields, 
e.g., considering or not s tellar rotation, could be ~0.2dex 
at solar metallicity (c.f., IRomano" et aI.ll2010T ). Thus, con- 
sidering quadratically the two effects, the absolute abun- 
dance level could vary by up to ~0.6dex without invoking 
infall/outflow processes. In any case, the approach presented 
here should be considered as a first step toward better under- 
standing the interplay of the three observables; future work 
will be aimed at comparing our current results with models 
having gas inflow and outflow and with different IMFs and 
stellar yields. 



4.2 The SFMS and the MZR 

Here we compare the model results with the observed SFMS 
(sSFR vs Af s tar relation) and the MZR (mass-metallicity 
relation). We also show the models superimposed on the 
Fundamental Plane identified in Paper I. To better visu- 
alize the interplay among the observables, and isolate dif- 
ferent regions of parameter space, we have divided the 
models and the data into two regimes of SFR: the "pas- 
sive" one with SFR<0.6Mq yr -1 , and the "active" one with 
SFR^O.6M yr" 1 . This threshold was defined in Paper I on 
the basis of the onset of deviations from the SFMS and the 
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Figure 3. Specific SFR vs. stellar mass divided into four ranges of SFR: SFR < 0.06 Mgyr -1 , 0.06 < SFR < O.6M yr-\ 0.6 < 
SFR < 6 Mq yr" 1 , SFR > 6 Mq yr" 1 . LVL (11HUGS) galax ies are shown as small open or filled circles (filled when SFR^0.6 M Q yr" 1 ), 
with different colors corresponding to Hubble type T as in lLee et al ^ 8 cyan, 5 ^ T < 8 blue, 3 ^ T < 5 magenta, and 

T < 3 red. Red X corresponds to the dlrr sample; BCDs are shown as (blue) squares (Fumagalli+Hunt samples) and (purple) triangles 
(Engelbracht); Green Peas are given by small (green) filled circles; LCGs by (green) +; LBGs at z ~ 1 (Shapley), z ~ 2 (Erb, Shapley) 
are shown as 6-pronged asterisks, and at z ~ 3 (Maiolino, Mannucci) as 8-pronged asterisks. 6-sided open stars show the two dwarf 
"prototypes", SBS 0335— 052 and IZwl8. Solid symbols show those BCDs with SFR^0.6 Mq yr" 1 . Models are coded by line type and 
color as in Fig. [2] The corresponding colour codes for the stellar masses (SE/RE, respectively) are: green/dark violet for 10 7 , 10 s Mq, 
light bluc/magenta for 10 9 , 10 10 Mq, and blue/red for 10 11 Mq. The evolution times are marked along the tracks: an inclined tick 
indicates 100 Myr (only for the RE galaxies), and a horizontal tick 1 Gyr. Time increases downward and to the right. The large + sign is 
a fiducial point (sSFR = 10~ 9 yr -1 , A/star = 10 9 Mq) which guides the eye to the differences in the main loci of data points and models 
in the different panels. 
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MZR; for our sample, these deviations become significant at 
all redshifts for SFR>0.6 Mq yr _1 , roughly the equivalent 
of NGC4826 (M64, the Black-Eye Galaxy), 1.5 times M33 
or the "dwarf starburst" NGC 3077^01 twice the SFR of 



M31 JPale et al.ll2012l;IVerlev et al.ll2007l ; iMeier et al.ll200ll ; 

iTabatabaei fc Berkhuiisenll2010l 'l. We have further divided 
these two regimes into two sectors, in order to better sep- 
arate and identify the progression from passive to active 
modes. In all of the following plots, we show models at times 
from 1 Myr to a Hubble time, as given in Fig. [2] 

4-2.1 The main sequence of star formation 

Figure [3] shows the SFMS with sSFR plotted against stellar 
mass. Galaxies characterized by a passive SF mode appear 
in the upper panels, while active galaxies are in the bot- 
tom ones. Models of RE and SE galaxies, with their mass 
"families" and three different MC density regimes, are also 
shown and divided into different SFR regimes. Evolution 
timescales are marked along the tracks: an inclined tick in- 
dicates 100 Myr (only for the RE galaxies), and a horizontal 
tick 1 Gyr. The + marks an arbitrarily defined fiducial po- 
sition in order to guide the eye to the displacements in the 
various panels as a function of SFR (see below). 

The least massive most quiescent galaxies (see upper 
left panel) have very low sSFRs, and are modeled only by 
the very end points of SE galaxies. As SFR increases (upper 
right), galaxies are more massive and the data are well re- 
produced by SE models with either low masses and compact 
MCs (the "intermediate" SF mode), or higher masses and 
diffuse clouds (the passive mode). Data in the bottom left 
panel (the onset of the active SF mode) , are well modeled by 
intermediate-mass SE or RE galaxies with dense or hyper- 
dense clouds, but also by the most massive SE galaxies (or 
RE galaxies) with diffuse clouds. As discussed in Sect. [3731 
both galaxy mass and cloud density contribute to SFR and 
evolutionary timescales. 

The most extreme starbursts in the bottom right panel 
are reproduced only by RE galaxies and the most massive SE 
galaxies with hyper-dense and compact clouds. Moreover, 
such high SFRs are achieved by the most massive RE models 
and the most massive hyper-dense SE only up to an age 
of ^ 200-500 Myr from the onset of the SF episode (see 
Fig. [2}. Thus galaxies with these high SFRs are expected 
to be "young", i.e., "captured" in a relatively early phase 
of their evolution. However, the RE models have assembled 
almost all their stellar mass (and metals) already at these 
early times, but the SE models will continue to build up 
their stellar and metal content for some time to come. 

The + in Fig. [3] represents a fiducial point (arbitrarily 
defined), with the aim of guiding the eye to the changes in 
loci among the panels. In the upper left panel, the galaxies 
fall below the fiducial by almost an order of magnitude in 
both axes; in the upper right and lower left, the galaxies 
(and models) are roughly coincident; in the last panel both 
galaxies and models fall above the fiducial by a factor of 5 
or so in both sSFR and stellar mass. Many LCGs and Green 
Peas at z ~ 0.3, and some BCDs at z ~ 0, are coincident 
with the LBGs in the AMAZE and LSD samples at z ~ 3. 
These galaxies are modeled by systems in early stages of 
their evolution, and with high SFR. Hence, at least in this 
parameter space of SFR and M sta r, there are local galaxy 



populations which mimic the properties of galaxies at high 
redshift. They are distinguished by their relatively high SFR, 
perhaps resulting from a similar selection effect for both the 
low- and high-redshift samples. 



4-2.2 The mass-metallicity relation 

The MZR is shown in Fig. [4] divided into four SFR regimes 
as in Fig. Low-mass galaxies with relatively high metal- 
licities (top left panel) are not well reproduced with our 
MICE models. This is probably because our SE models re- 
quire more than a Hubble time to complete the stellar mass 
assembly and metal production with the available fuel. On- 
going accretion could help here, but it would also be nec- 
essary to shorten timescales in order to accelerate the gas 
consumption. As SFR increases (top right), intermediate- 
mass SE galaxies are able to well reproduce the data, either 
more massive models with diffuse clouds or less massive ones 
with compact or hyper-dense clouds. 

Galaxies forming stars in the active mode (bottom left 
panel) are well reproduced by the more massive SE galax- 
ies and the less massive RE galaxies. The most massive SE 
galaxy (10 1 Mq) with diffuse clouds is consistent with the 
most metal-rich observations, and also the most metal-poor 
ones. Interestingly, the latter appear to lie along the same 
evolutionary path as the former and may evolve into them. 
However, an active SF mode with compact or hyper-dense 
clouds is necessary to reproduce the observations at lower 
masses: either SE galaxies for masses 10 s - 10 9 M Q , or RE 
galaxies for lower masses. 

Extreme starbursts (bottom right panel) require rapidly 
evolving models at young ages for such high SFRs; galax- 
ies must be younger than < 200-300 Myr (see Fig. [2J. The 
main population of the galaxies in this panel can be repro- 
duced only with the early phases of the evolutionary tracks 
of RE (starburst) galaxies. Only a single SE galaxy model, 
the most massive one with hyper-dense molecular clouds, 
appears in this panel. Below, we explore how the gas con- 
tent of the models changes with SFR (and age) in order to 
better understand the evolutionary phase of these galaxies. 



4.3 The fundamental plane of metallicity, SFR, 
and stellar mass 

Our models have been rather successful at reproducing the 
observed trends in the SFMS and MZR. However, in Paper I, 
we introduced a Fundamental Plane (FP) which enabled an 
accurate (0.17dex) determination of O/H given SFR and 
Mstar- Here we show that our models also well reproduce 
the FP, and occupy the same spread in parameter space as 
the observations. 

Figure [5] shows the Fundamental Plane that emerged 
from the Principle Component Analysis (PCA) of Paper I. 
The two upper panels show the edge-on views of the plane, 
and the bottom panel the face-on one. Our models are su- 
perimposed on the data and are consistent with the loci of 
the data in all three projections. In all three plots, time in- 
creases downward. 
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Figure 4. Nebular oxygen abundance vs. stellar mass M a t a r divided into regions of SFR. The four regimes we consider are: SFR ^ 
0.06M Q yr~ 1 , 0.06 < SFR < O.6M yr- 1 , 0.6 < SFR ^ 6M yr- 1 , SFR ^ eAfeyr" 1 . Symbols are as in Fig. [3] Models are coded by 
line type and color as in Fig. [3] diffuse as dotted lines, compact as dashed, and hypcr-dense as solid lines. The corresponding colour 
codes for the stellar masses (SE/RE, respectively) are: green/dark violet for 10 7 , I0 V M . . light blue/magenta for 10 9 , 10 10 Mq, and 
blue/red for 10 11 Mq. The evolution times are marked along the tracks: an inclined tick indicates 100 Myr (only for RE galaxies), and 
a horizontal tick 1 Gyr. Time increases upward and to the right. The large + corresponds to a "fiducial" value of iWstar = 10 9 Mq, and 
12+log(0/H) =8.0, intended to guide the eye to differences among the panels as in Fig. [3] 



5 GAS FRACTIONS AND METAL 
ABUNDANCE 

Most previous interpretations of the shape of the MZR 
have relied on gas inflow and outflow, with the assumptions 
that outflowing gas is chemically enr i ched and that infalling 
gas is not (e.g.. iTolstov et all [20091 : iMannucci et al.ll20ld ; 



iDave et alj|2012l : lYates et all [20121 : iDaval et ail I20121 ). As 

already discussed, our MICE models are "closed-box", and 
thus neglect the effects of gas inflow and outflow. Neverthe- 
less, as we have shown in previous sections, the predictions of 
our models are in good agreement with the SFMS and the 
MZR through the active/passive dichotomy of SF modes 
and the resulting evolutionary timescales. Here we exam- 




Figure 5. Different projections using the 3 PCs found by the PCA. Galaxies are coded as shown in the lower right (empty) panel. 
The top left and right panels show the orthogonal "edge-on" views of the plane; the bottom panel shows the plane face-on. x\ = 
12+log(0/H)-( 12+log(0/H)), x 2 = log(SFR)-{log(SFR)) (Mq yr^ 1 ), and xz = log(M atar )-(log(M star )) (M ). (12+ log(0/H)> = 
8.063; (log(5Ffi)}= -0.594; (log(M st ar)> = 8.476. PCI = 0.12 x x + 0.75 x 2 + 0.65 x 3 ; PC2 = -0.31 xi - 0.65 x 2 - 0.69 x 3 ; PC3 = 
-0.94xi - 0.11x2 + 0.31x3. The axes arc expanded in the top panels to exaggerate the variations in PC3, relative to the much larger 
dynamic ranges in PCI and PC2. MICE models are superimposed similarly to previous figures, but here there are no temporal tick 
marks. In all three panels, time increases downward. 



ine the MICE model predictions for the molecular gas mass 
fractions relative to stellar mass, and compare them with 
observations of local and high- 2; samples with measured H2 
and stellar masses. 

The time evolution of the H2 content ("cloud" mass) 
and the cloud-to-stellar mass ratio is shown in Fig. [6] The 
right panel of this figure, as in Fig. [2] shows the degeneracy 
between the track of the SE galaxy with M ga j = 10 11 Mq 
and that of the RE galaxy with M ga i = 10 8 Mq. The scaling 
factors conspire to make these two galaxies behave similarly 
in any quantity which is normalized (i.e., O/H, H^/Mjtar)- 
However, as before, the different timescales for the different 
SF modes make the gas fraction relative to stars (and chemi- 
cal abundance) non-unique tracers of evolution. A low-mass 
RE spheroid can achieve similar metallicity and gas mass 
fraction as a more slowly evolving, but more massive, SE 



galaxy. The gas fraction of the rapidly evolving spheroid 
will be higher than that of the more slowly evolving galaxy, 
as long as we observe the former at an earlier evolutionary 
phase than the latter. 

In Fig. [Jj we have plotted the ratio of molecular cloud 
to stellar mass vs. Af s t ar , separated into four panels of 
SFR as in previous figures. Also shown are s everal sam- 
ples including the BCDs from iFumagalli et al.l |2010h (also 
considered i n the previous analyses), spiral galaxies from 
iLerov et ail (|2008|). galaxie s from the COLDGASS survey 
by Saintonge et al.l ( 2011alfl), and sub-m i llimet re galaxies 
(SMGs) at z ~ 2 - 3 from lHainline et al l (|201ll r1. Conver- 
sion factors Xco to derive H2 mass ( "cloud mass" in Fig. [7} 



4 There is only one sample in common with the SFMS and MZR 
analysis, because of the current unavailability of molecular gas 
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Figure 6. Predictions of H2 mass (left panel) and ratio of H2 to M star (right panel) vs. time in the MICE models. As in Fig. [2] we show 
only models of two representative galaxies, 10 s M@ (green/dark violet) and 10 11 Mq (blue/red), for the SE/RE types, respectively. The 
diffuse (dotted lines), compact (dashed), and hyper-dense cases (solid lines) are coded by line type. In these models, the cloud-formation 

1/3 

efficiency and SFR efficiency both scale with cloud mass as: Af MC . 



are given in the respective papers, and vary with galaxy 
type. 

The main problem with comparing molecular content 
of metal-poor starbursts (lower left and lower right pan- 
els in Fig. [Jj lies with the intrinsic difficulty of detecting 
molecules in the c o ol gas phase in a metal-poor ISM (e 
iTavlor et all Il998l ; iBarone et all l200fj| ; iLerov et al] |200~ 
Moreover, the conversion of CO intens ities to H2 mass at 
low metallicities i s very uncertain (e.g. , IBolatto et al ]|2008l : 
ILerov et al] 120091; IWolfire et all 120101 ; ILerov et al] l201ll ). 
While IBolatto et all ( 20081 ) find Xco not to vary signifi- 
cantly with metal abundance in Local Group galaxies, at 
metallicities 12+log(0/H) ,$,8.2 t here is evidence for Xc o 
factors of 10 to 20 times higher (jLerov et all 120091. l201lT). 
For t he BCDs with CO observations in iFumagalli et al] 
(2010). we have used a conse rvative value, Xco = 4x 
10 20 cm" 2 [Kkms" 1 ]" 1 , (e.g.. IBolatto et all 120081 ); hence, 
the estimates of H2 mass in the metal-poor BCDs shown in 
Fig. [Jj could be highly underestimated (and in fact are out- 
liers because of their low inferred cloud mass/stellar mass 
ratios). 

The MICE model predictions of the ratio of molec- 
ular gas mass to stellar mass shown in Fig. [Jj are in 
fairly good agreement with the observations, excepting the 
COLDGASS upper limits in the upper left panel with 
SFR<0.06 Mq yr -1 . This consistency is perhaps surprising, 
given that our models are closed box, but strongly suggests 
that other parameters, such as dynamical timescales result- 
ing from compact size and high density, play an important, 
if not primary, role in governing gas fractions. 

Comparison with Figs. [3] and [4] indicate that the mod- 
els that reproduce the scaling relations between SFR, stellar 



observations in galaxies with measured stellar masses, SFR, and 
metallicities. 



mass, and O/H in the upper panels are relatively gas poor. 
The upper left panel in Fig. [3] is not well reproduced by our 
models because the galaxies tend to be more massive than 
expected from the sSFR; this is perhaps propagating into 
Fig. [Jj where galaxies are again too massive than would be 
predicted from their low SFR and low cloud mass fraction. 
The upper right panel is apparently occupied by either RE 
RE galaxies or more slowly evolving galaxies which, in any 
case, are at the end of their lifetimes, having already con- 
verted most of their gas into stars. As the SFR increases 
(lower left panel), gas content also increases, and for the 
most extreme case (lower right panel), the mass ratio of 
molecular clouds to stars is ~1, during the time when the 
galaxies also satisfy the scaling relations between metallic- 
ity, stellar mass, and SFR. The data for SMGs are entirely 
consistent with the models, and in fact have observed gas-to- 
total mass fractions [clouds/ (stars+clouds)] of 50% or more 
ijTacconi et al.ll2006l ). 

With our MICE models we are able to link the gas 
content and the SFR with the evolutionary timescales. The 
galaxies we call active are observed close to the peak of the 
SFR with a still high gas content, and are relatively chem- 
ically unevolved. During that phase of their evolution they 
are still gas-rich star-forming galaxies but will rapidly move 
toward the end of their SF phase. In fact, their evolution 
can be reproduced using chemical evolution models for RE 
(starburst) galaxies. On the other hand, passive galaxies are 
evolving at a very slow rate, so that their gas content and 
SFR will be relatively low during most of their lifetime. 

Recently, IZahid et al] (|2012l ) presented a sample of 
metal-rich low-mass galaxies, a galaxy population comple- 
mentary to the metal-poor more massive starbursts dis- 
cussed here. They interpreted the relatively high metallic- 
ities in their sample as due to the reduced gas content of 
the galaxies, such that the more metal-rich galaxies at a 
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Figure 7. Ratio of molecular cloud mass and M s tar vs. stellar mass M s tar divided into regions of SFR. The four regimes we consider are: 
SFR ^ 0.06 M fl yr-' , 0.06 < SFR sC 0.6 Mq yr" 1 , 0.6 < SFR ^ 6M Q yr " 1 , SFR ^ 6M ff y r-'. Squares (blue) correspond to the BCDs 
fromlFumagalli et al.l l|2010l ). triangles (green) to the spiral galaxies from lLerov et al and 6-pronged asterisks to the SMGs fr om 

lHainline et al. I j201lh . Turquoise + and down-arrows indicate the COLDGASS molecular gas content from ISaintonge et al. | J2011a| [bh. 
Models are coded by line type and color as in Fig. \3\ The evolution times are marked along the tracks: an inclined tick indicates 100 Myr 
(only for RE galaxies), and a horizontal tick lGyr. Time increases downward and to the right. The large + corresponds to a "fiducial" 
value of cloud mass/M s tar of 0.1 and M s t a r = 10 9 Mq, intended to guide the eye to differences among the panels. 



given stellar mass are more gas poor. In fac t , the metal- 
rich low-mass galaxies studied bv lZahid etafl l|2012h would 
occupy the upper panels in Fig. [7| The interpretation of 
IZahid et al.l (|2012l ) is that at a given stellar mass, metal-rich 
galaxies have low gas fractions, and that the low gas frac- 
tions are largely responsible for the high gas-phase oxygen 
abundance observed. Instead, our models suggest that the 



highest metallicities at a given mass are not a consequence of 
the gas fraction, but that both gas fraction and metal abun- 
dance are consequences of evolutionary time scales, driven 
by the initial conditions in passive and active galaxies. Thus, 
our models also appear to ver ify the empirical conclusion 
reached by IZahid et all (|2012l ) that such galaxies are gas 
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poor relative to their metal-poor counterparts at the same 
stellar mass. 



6 DISCUSSION 

A significant fraction of galaxies in our sample deviates 
from the MZR and SFMS in similar ways, independently 
of redshift. These galaxies are apparently too massive for 
their metallicity, or alternatively, for their mass they are 
too metal-poor and have a SFR that is too high. The com- 
mon characteristic of these galaxies, the low-metallicity star- 
bursts, is the way they are selected from larger samples. 
Locally, metal-poor starbursts tend to come from objective- 
prism, emission line surveys; they are objects with strong 
H/3, [Oiii]A5007, or Ha lines, and thus moderate dust extinc- 
tion and high SFRs. At z 1, Lyman break selections, based 
on rest-frame UV colors, pick out star-forming galaxies with 
similar properties, namely relatively low dust extinction and 
relatively high SFRs fe.g.. lShapley||20riT ). Such galaxies also 
tend to be dominated by the most recent burst of star forma- 
tion, with times from burst of a few tens of M yr to 1 Gyr 
|Shaplev et al. I l2005bl ; iHaberzettl etalll2012h . At low red- 
shift, this is a result of the selection of very high emission- 
line equivalent widths either in Ha or in [Om] which implies 
that ionized gas makes a significant contribution to the to- 
tal light (implying a young age). At high redshift, this is 
a function of the color selection for strongly star-forming 
galaxies that can be finely tuned to select galaxy popula- 
tions at different redshifts. Hence, the low- and high-redshift 
LMSs, which deviate from the SFMS and the MZR in simi- 
lar ways, have similar properties because of the way they are 
selected. In fact, together they form a Fundamental Plane of 
O/H, SFR, and M sta r, which is apparently independent of 
redshift. 

The galaxies which deviate from the SFMS and MZR 
are well modeled by the active mode of star formation. The 
starbursts in the lower right panels of Figs. [3] and 2] are ob- 
served during the extreme starburst phase of their evolution 
which is rapid and intense, lasting up to £ 200-500 Myr. On 
the other hand, the galaxies which follow the general trends 
in the scaling relations form stars in the passive mode. The 
models suggest that they are observed during a more ad- 
vanced long-lasting phase of their evolution, which endures 
most of their lifetime. 

What are the drivers of the active and passive modes of 
star formation? There are two related aspects of our models 
which differentiate active and passive SF modes: their in- 
trinsic characteristics in terms of size (RE vs. SE) and/or 
density of their star forming regions (diffuse, compact, and 
hyper-dense clouds), and their evolutionary age. The evolu- 
tionary age indicates the time from the last main episode of 
SF: for quiescent galaxies which have had an initial single 
main episode of SF, this might be comparable to a Hubble 
time, while for active galaxies it might be the time elapsed 
after the last major episode of SF. Thus when we refer to 
active galaxies as young objects, it means that they had 
their major episode of SF relatively recently, up to several 
hundreds of Myr ago. 

The predictions of our models are in agreement with re- 
cent analy ses of the structure of galaxies up to z ~ 1 along 
the SFMS jWuvts et al.ll201ll '). At all stellar masses and red- 



shifts, they find that the largest disk-like galaxies lie along 
the SFMS; these would be the galaxies with the longest dy- 
namical times. The star-bursting outliers are more centrally 
concentrated (compact) than the galaxies along the main se- 
quence; these would have shorter dynamical times and thus 
be more rapidly evolving. Thus, our models correctly pre- 
dict that compact RE active models reproduce the starburst 
outliers of the SFMS, and that SE galaxies, with their larger 
sizes and slower evolution, lie along it. 

An active SF mode is more readily associated with 
mergers and interactions than a passive one; in fact, the ini- 
tial conditions we attribute to active star formation may be 
induced by mergers and tidal compression. However, merg- 
ers are not the only way to foster a highly star- forming mas- 
sive galaxy at z ~ 2 (e.g., iGenzel et al.|[200r} ). In fact, we 
hypothesize that the initial conditions for the active and 
passive modes of SF are governed by stochastic processes, 
which can occur even in isolation. 

Although our MICE models reproduce the scaling rela- 
tions of metallicity, stellar mass, and SFR rather nicely, they 
have limitations; as already discussed, the main one is they 
do not take into account either inflowing or outflowing gas. 
Other attempts to explain the MZR (and FMR) have been 
based on the transfer of gas, i.e., accretio n of metal-poor 
gas and outflow of metal-enriched mater i al (iMannucci et alj 
l2010l : lYates et alj|2012l ; iDave et alj|2012l : iDaval et alj|2012h . 
R ecent theore t ical w ork, such as the "equilibrium approach" 
of lDave et ail (|2012h . consider stochastic rates of inflow and 
outflow that regulate SFR and the accumulation of stellar 
mass and metallicity. Such models are quite successful in 
approximating the shape of the MZR and its dependence on 
SFR, at least down to 12+ 1or(0/H)~8, SFR~1 M© yr -1 , 
and AW ~ 10 s ' 5 M© (|Yates et alJl2012T ). It is not clear, 
though, whether they are consistent with the lower metallic- 
ities and higher SFRs in the deviant metal-poor starbursts. 

Another limitation of our models is the way we approx- 
imate the onset of the molecular phase; the atomic gas must 
become molecular and form clouds before it can fuel star 
formation. Our models assume that molecular clouds form 
at a constant efficiency, with a dependence on atomic gas 
mass, and are disrupted by cloud-cloud collisions. More com- 
plex models take into account the self-shield ing effect of dust 
on the UV interstellar radiation field (e.g.. iKrumholz et alj 
120091 : iKrumholz fc Gnedidl201lh . Although our models are 
successful in predicting the trends of the SFMS and the 
MZR, and observed gas fractions as in Fig. [7] they may be 
oversimplified. A more sophisticated treatment of the tran- 
sition to molecular gas will be dealt with in a future paper 
(Schneider et al., in preparation). 

Despite the limitations of the MICE models, they suc- 
cessfully predict the observed trends in scaling relations over 
a wide range of O/H, SFR, and Af star . Their predictions are 
also in good agreement with observed molecular gas mass 
fractions and SFR. It is likely that some of the processes 
driving galaxy evolution are captured by our models, namely 
the formation of a dominant molecular phase, together with 
fast or slow evolutionary timescales depending on stochastic 
initial conditions. It is also likely that the physics behind 
galaxy evolution requires some gas exchange with the envi- 
ronment, and a regulation mechanism t hat drives an equ i- 
librium between inflow and outflow (e.g.. |Pave et alj|2012h . 
A combination of key features of both classes of models may 
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be necessary to understand the mechanisms which drive evo- 
lution of galaxies in the early universe. 



7 CONCLUSIONS 

From a sample of ~ 1100 galaxies spanning a wide range 
of O/H, SFR, and M sta r, we have defined a class of low- 
metallicity starbursts which deviate from the SFMS and the 
MZR in the same way, independently of redshift. The sam- 
ple taken as a whole (although excluding the galaxies with 
Mrtar^ 3x 10 10 Mq) defines a Fundamental Plane whose ori- 
entation is defined primarily by SFR and Af s t a r, and whose 
thickness is governed by O/H. Our models show that this 
plane is populated by galaxies and by models with different 
modes of star formation: an active starburst mode that is 
more common at high redshift, and a more passive, quies- 
cent mode that predominates in the local universe and that 
defines the main trends in the scaling relations. Summariz- 
ing the results shown in Figs. [3l [4] and [7] galaxies follow- 
ing the main trends in the SFMS and MZR form stars in 
the passive mode. They are well modeled by slowly evolving 
galaxies with only moderately dense molecular clouds (ei- 
ther "diffuse" or "compact" MCs) and a relatively low gas 
content. These galaxies spend most of their life in this qui- 
escent, passive, evolutionary phase reaching asymptotically 
the conversion of their gas into stars in more than a Hubble 
time. We find that observations are consistent with the end- 
phases of the evolution of these galaxies predicted by our 
models. On the other hand, the galaxies which deviate from 
the SFMS and MZR form stars in the active mode, typical 
of gas-rich RE galaxies with masses from 10 8 to 10 11 Mq, 
or gas-rich SE galaxies with compact or hyper-dense clouds 
with masses of ~ 10 11 Mq. In the extreme cases of very 
high SFR, the evolutionary phase is short, suggesting that 
we are observing these galaxies just after the start of their 
most recent SF episode. 

The existence of a Fundamental Plane and the success of 
our modeling of active/passive star formation implies that 
the SFMS and MZR do not really change with lookback 
time, but rather that the galaxy populations that define 
them change with redshift. The "evolution" in the scaling re- 
lations is a result of selecting those galaxies which are most 
common at a particular redshift. Equivalently, the scaling 
relations and the deviations from them result from different 
modes of star formation: an active starburst mode that is 
more common in the distant universe and a more passive 
quiescent one that dominates nearby. 
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